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Using the domino dynamo model, we show how variations of the heat flux at the core-mantle 
boundary change the frequency of geomagnetic field reversals. In fact, we are able to demonstrate 
the effect known from the modern 3D planetary dynamo models using an ensemble of interacting 
spins, which obey equations of the Langevin type with a random force. We also consider applications 
to the giant planets and offer explanations of some specific episodes of the geomagnetic field in the 
past. 
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Generation of the planetary magnetic fields is a subject 
of the dynamo theory, which describes successive trans- 
formations of thermal and gravitational energy, concern- 
ing compositional convection, to energy of kinetic mo- 
tions of the conductive liquid and then to the energy 
of the magnetic field [l]. Modern dynamo models in- 
clude partial differential equations of thermal and com- 
positional convection as well as the induction equation for 
the magnetic field, which some reasons should be three 
dimensional [2]. Although, due to the finite conductivity 
of the Earth's mantle, observations of the geomagnetic 
field at the Earth's surface are bounded by the first thir- 
teen harmonics in the spherical function decomposition, 
one needs small-scale resolution down to ^ 10~^ L, to 
provide the necessary force balance in the core. Here 
L ~ 2.2 lO^m is the scale of the liquid core. This dif- 
ficulty is caused by the huge hydrodynamic Reynolds 
number Re ~ 10^ as well as by the strong anisotropy 
of convection [3] due to the geostropic state in the core 
Convection in the core is cyclonic. The cyclones and 
anticyclones are aligned with the axis of rotation, and 
their scale is much smaller than their length. As a result 
one needs very efficient computer resources to produce 
regimes in the desired asymptotic limit required for geo- 
dynamo simulations with grids 128^ and more, which is 
a challenge even for modern supercomputers. 

In spite of these technical problems, 3D dynamo mod- 
els successfully mimic various features of the modern and 
ancient magnetic field including the reversals One of 
the important results of the dynamo theory is that the 
frequency of the reversals depends on the spatial distri- 
bution of the heat fiux at the outer boundary of the liq- 
uid core jif. In particular, the authors have shown that 
the increase of the heat fiux along the axis of rotation 
leads to increase of the axial symmetry of the whole sys- 
tem and stops reversals. In a sense, the thermal trap of 



reversals occurs. On the contrary, the decrease of the 
thermal fiux at high latitudes leads to the chaotic be- 
haviour of the magnetic dipole accompanied by frequent 
reversals, which is closely connected with the upsetting of 
the geostrophic balance and predominance of the radial 
(in an incompressible medium in which the parameters 
depend only on the radius) Archemedean forces. It looks 
attractive to obtain this result using toy dynamo models 
(such as the Rikitake and Lorentz models), which can 
provide extensive statistics and obviousness of the re- 
sults, and simulate practically instantly, using just home 
PC, the number of reversals and excursions of the same 
order as known from paleomagnetism. The random force, 
which imitates the small-scale unresolved fiuctuations, 
was included in some of these models 0. 

For this reason we have selected the domino dynamo 
model [8, 9], which is an extension of the Ising-Heisenberg 
XY-models of interacting magnetic spins, for more de- 
tails of the history of the problem and classification refer 
to 

The main idea of the domino model is to consider a 
system of N interacting spins Si, i = 1 . . .N, in media 
rotating with angular velocity ft = (0, 1) . The spins 
are located over a ring, are of unit length and can vary 
angle 9 from the axis of rotation on time t, so that Si = 
(sin^i, cos^i). Each spin Si is forced by a random force, 
effective friction, as well as by the closest neighbouring 
spins Si_i and S^+i. 

Following [8], we introduce kinetic K and potential U 
energies of the system: 
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7 5^(12-S.)' + A5^(Si-Si+i). 

i=l n=l 
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The Lagrangian of the system then takes the form £ = 
K — U. Making the transition to the Lagrange equations, 
adding friction proportional to 9 and the random force 
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X, 

d dC ^dC _ ■ ex 

leads to the system of Langevin-type equations [11 : 

'Oi — 27 cos 6i sin + A |^ cos 9i ^ sin 6i-i + sin ^i+i 
sin^i^cos^i-i + cos^i+i^l + re^i + ^ = 0, 
^0 = ^iVi ^7V+i = ^1, i = 1 . . . iV, 

(3) 

where 7, A, e, r are constants. The measure of syn- 
chronization of the spins along the axis of rotation 
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choice of parameters) to paleomagnetic observations [Tl| , 
including the periods of the random and frequent rever- 
^2) sals. It was shown in [1] that large fluctuation of a single 
spin can be successively transferred to the neighbouring 
spins, which fluctuate until all spins reverse their polar- 
ity. This dynamo effect was an inspiration for the name 
"domino mode" .1 

In all simulations we have used, similarly to [8], val- 
ues of parameters 7 = —1, A = —2, n = 0.1, e = 0.65, 
T = 10~^, iV = 8, with random normal Xii with zero 
mean values and unit dispersions, so that Xi was updated 
at every time step equal to r. As follows from [8jthe evo- 
lution of M depends slightly on the form of the random 
forcing x ^-^^d the remaining parameters can be easily 
selected in such a way as to provide similarity with ob- 
(4) servations. The typical behaviour of the axial dipole is 



will be considered to be the total axial magnetic moment. 
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FIG. 1: Evolution in time of M: (a) forthe purely magnetic 
system {C^ = 0), (b) = 0.5, (c) = -0.5 for ip = 
- cos^ 6, (d) = 10, (e) = -9 for V = - cos^ 29. At (d) 
the thick line corresponds to M and the thin line to Me- 



Integrating (jSj in time for small N one can arrive at 
quite diverse dynamics of M, very close (for some special 
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FIG. 2: Plots of cos Oi with i = 1 . . . N , versus M for reversal 
in Fig. la during the time interval t = 10^-3 • 10"^. 



presented in Fig{T^, where we observe several reversals at 
irregular time intervals. This process is accompanied by 
a short drop of M to nearly zero and followed by rapid 
recover, which in geomagnetism is referred to as excur- 
sions of the magnetic field. One can find a thorough 
analysis of this system in jsf ; here we only emphasize one 
important in geomagnetism point. Up to now observa- 
tions do not clearly indicate, if the geomagnetic dipole 
just rotates during the reversal without decreasing the 
amplitude, orif it decreases and then recovers with the 
opposite sign. To check these scenarios we present the 
evolution of the individual spins during the reversal, see 
FiglU There are quite large deviations of the individual 
spins from the mean value M at the moment of reversal 
(M = 0). This means that the decrease of M is caused 
by the desynchronization of the spins rather than by the 
coherent rotation of the spins. Note that the typical time 
of the reversal is much larger than the time step. The 
points on the lines correspond to every lO*'^ time step in 
the simulations. The other point is that the minimal time 
^ N t/2 of propagation of the disturbance from spin Si 
through half the circle is also smaller than the typical 
reversal time. This scenario is also supported by the 3D 
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simulations, where the spots of the magnetic field with 
opposite polarities co-exist at the core-mantle boundary 
during the reversals. 

We now extend the concept of the spin from the purely 
magnetic system to the whole cyclone system, including 
its hydrodynamics, and we introduce correction ^ to po- 
tential energy U ([1]), which takes into account the het- 
erogeneity of the thermal flux. For the purely magnetic 
problem, it corresponds to the electromagnetic interac- 
tion of the cyclone with the conducting mantle, the D" 

layer. The new effective force Fi = — — -- , whose infiu- 

ence on the behaviour of M{t) will be considered in the 
rest of the paper, then appears in the right-hand side of 
®. 

Let ^{t, 6) = t/j{t, 9), where is a constant, and 
the spatial distribution of the potential is given hy ip = 
— cos^ 6. Then, > corresponds to the stable state 
in the polar regions, ^ = 0, tt, and the appearing force 
F = — sin 29, acting on the cyclones, is directed towards 
the poles. This regime corresponds to the increase of the 
thermal flux near the poles that causes the stretching 
of the cyclone along the axis of rotation. In Fig{T] we 
demonstrate the effective influence of F on M. 

The increase of the thermal flux along the axis of rota- 
tion leads to the partial suppression of the reversals of the 
field, Fig{T]3. Note that, for the chosen potential barrier 
tjj, the dependence of F is equal to the 7-term in (jS)): the 
increase of the thermal flux at the poles leads to the effec- 
tive increase of rotation and amplification of geostrophy, 
caused by the rapid daily rotation of the planet. Our re- 
sults are in agreement with the 3D simulations, see Fig. Id 
in . The further increase of [C^ = 2) leads to the 
total stop of the reversals. It is more interesting that, 
using even larger > 10, one arrive at regimes with 
a nearly constant in time |M| < 1 defined by the initial 
distribution of Sj. In other words, the super flux at the 
poles can fix the spins which are still not coherent. There 
is some evidence [l2j, |l3| that the geomagnetic dipole in 
the past couldhave migrated from the usual position near 
the geographic poles to some stable state in the middle 
latitudes. Within the framework of our model, we can ex- 
plain this phenomenon by the thermal super flux at the 
poles. Later we will discuss some other scenarios which 
yield similar results. 

For negative C^, when the geostrophy breaks due to 
the relative intensification of convection in the equatorial 
plane, we get the opposite result, see Fig(TJ:;: the regime 
of the frequent reversals observed in Fig.lc in In 
this case force F is directed from the poles and the equi- 
librium point at the poles becomes unstable. The new 



minimum of the potential energy at the equator leads to 
the appearance of a new attractor, so that for = — 5 
the axial dipole fluctuates with zero mean value and max- 
imal amplitude M ^ 0.4. This state corresponds to the 
equatorial dipole 

1 ^ 

Me{t) = -J2^'''^^(t), (5) 

i=l 

so that I Me I ^ 1 and does not undergo reversal. Similar 
behaviour of the magnetic dipole is observed on Neptune 
and Uranus; for more details see, e.g. [l^ . 

Now we consider tp = — cos^ 29 for which the cor- 
responding force F = —2 sin 4^ changes sign in each 
hemisphere, see the example with the second-order zonal 
spherical harmonic ^ P2, where Pf^ is associated Legen- 
dre polynomials, in Fig.le in [6[. 

For > 0, the potential barrier is in the middle lati- 
tudes, which prohibits the reversals as observed in Fig.le 
in Pi . In one of our runs only one reversal was observed 
at Ci/, = 2. The other important point is the existence 
of the stable point at the equator, where M = 0. We 
could then assume a regime with two attractors: near the 
poles and at the equator. This regime is really observed 
in FigllJl. The inverse transition was not observed. In- 
spection of the state with small M at the beginning of 
the run leads to a quite large estimate of the equatorial 
dipole amplitude Me, see Fig{TJl. Thus, in principle, this 
regime can also be related to the giant planets' dynamo 
as. 

The last example corresponds to the case when < 0. 
Here in addition to the attractors at the poles (related 
to rotation O), two new attractors appear at middle lat- 
itudes. The variation of C^, leads to regimes = — 1 
with frequent reversals, observed in Fig. If in More- 
over, we can get regimes, see Fig{T^, when the magnetic 
pole stays at high latitudes, however, |M| / 1 (the par- 
tial synchronization of the spins). There are some jumps 
to the unstable state M = with the dipole at the equa- 
tor. Note that the decrease of the spatial scale of ip leads 
to the increase of its amplitude C^. 

Obviously, as the interaction between the cyclones is 
certainly more complex [15] than in ([1]) and the mecha- 
nism of magnetic field generation in cyclonic turbulence is 
outside the scope of the domino model, our results should 
be used with caution. However, it is remarkable that even 
such a simple model provides quite a wide range of effects, 
known from very sophisticated models and observations. 
This, of course, does not exclude other explanations and 
scenarios of the considered phenomena. 



[1] G. Rudiger, R. Hollerbach. The Magnetic Universe: Geo- [3] P. Hejda, M. Reshetnyak. Effects of anisotropy in the 

physical and Astrophysical Dynamo Theory. 2004. Wein- geostrophic turbulence. Phys. Earth Planet. Int. 2009. 

heim. 177, 152-160. 

[2] C. A. Jones, Annu. Rev. Fluid Mech. 43, 583 (2011) [4] Pedlosky J 1987 Geophysical Fluid Dynamics (NY, 



4 



Springer- Verlag) 
[5] U.R. Christensen, J. Fluid. Mech. 470, 115 (2002) 
[6] G.A. Glatzmaier, R. S. Coe, L. Hongre, P. H. Roberts, 

Nature. 401, 885 (1999) 
[7] P. Hoyng, Astron. Astrophys. 272, 321 (1993) 
[8] N. Mori, D. Schmitt, A. Ferriz-Mas, J. Wicht, H. Mouri, 
A. Nakamichi, M. Morikawa. A domino model for geo- 
magnetic field reversals. (2011). [arXiv:lll0 .5062vl. 
[9] A. Nakamichi, H. Mouri, D. Schmitt, A. Ferriz-Mas, J. 
Wicht and M. Morikawa, Submitted to Mon. Not. R. 
Astron. Soc. (2012). arXiv:1104.5093vl. 
[10] H. E. Stanley, Introduction to phase transitions and crit- 
ical phenomena. 1987. Oxford. 
[11] J. A. Jacobs, Reversals of the Earth's Magnetic Field. 



1994. Cambridge. 

[12] A. V. Shatsillo, A. N. Didenko and V. E. Pavlov. Two 
competing Paleomagnetic directions in the Late Vendian: 
New data for the SW Region of the Siberian Platform, 
Russ. J. Earth Sci. 7, 4 ES4002 (2005) 

[13] A. Abrajevitch, Rob Van der Voo. Incompatible Edi- 
acaran paleomagnetic directions suggest an equatorial 
geomagnetic dipole hypothesis. Earth and Planetary Sci- 
ence Letters. 293 164 (2010) 

[14] I. Cupal, P. Hejda, M. Reshetnyak, Planet. Space Sci. 50 
1117 (2002) 

[15] A. Kageyama and T. Sato, Phys. Rev. E5, 4, 4617 (1997) 



